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Inclined  pileup  of  edge  dislocations  near  the  crack  tip 


Guido  D.C.  Dhondt  and  A.  Cental  Eringen 
Princeton  University,  Princeton  NJ. 


Abstract. 


The  formation  of  a  plastic  zone,  in  the  form  of  a  linear  array  of  inclined  edge 
dislocations  in  front  of  a  mode  I  crack  in  plane  strain,  is  studied.  The  direction 
of  the  force  field  on  an  edge  dislocation  due  to  the  crack  load  allows  for  the 
determination  of  the  plastic  zone  as  a  function  of  the  inclination  angle.  The 
dislocation  distribution  along  the  plastic  zone  is  found  from  force  equitibritim 
considerations.  The  results  show  the  presence  of  a  dislocation  free  zone  (DFZ) 
and  a  shrinking  of  the  plastic  zone  size  with  greater  angle  of  inclination.  The  latter 
angle  does  not  seem  to  influence  the  DFZ  size.  The  way  in  which  we  determined 
the  stress  field  allows  us  to  study  the  influence  of  defects  on  the  plastic  zone 
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Introduction 


Electron  microscope  studies  have  revealed  that  crack  propagation  involves  a  com¬ 
bination  of  modes.*1’2-3-4*5)  In  thin  areas,  where  the  state  of  stress  is  predominantly 
that  of  plane  stress,  the  crack  propagates  as  a  mode  IH-crack  and  emits  screw 
dislocations  in  a  plane  coplanar  with  the  crack.  In  thicker  areas  the  crack  behaves 
as  a  mode  1-crack,  and  either  it  blunts  by  emitting  edge  dislocations  on  planes 
45°  to  90°  inclined  to  the  crack  plane  or  propagates  suddenly  in  a  zig-zag  fashion. 

The  distribution  of  dislocations  in  the  plastic  zone  coplanar  to  the  crack  was  first 
treated  theoretically  by  Bilby,  Cottrell  and  Swinden  (BCS).*6)  The  pile-up  of  screw 
dislocations  in  front  of  the  mode  Ill  crack  was  later  studied  by  Majumdar  and 
Bums/7)  and  Chang  and  Ohr.**>  These  latter  studies  included  the  presence  of  a 
dislocation-free  zone  (DFZ).  Chang  and  Ohr  also  studied  the  inclined  pile-up  of 
screw  dislocations,  without  DFZ*9). 

Experimental  observations  on  the  plastic  zones,  ahead  of  the  crack  tips,  in  real 
materials  are  substantially  different  in  that  the  plastic  zones  are  spread  on  two 
thin  leaves  attached  to  each  of  the  crack  tips.  Here  we  idealize  this  by  conadrring 
one  or  two  symmetric  lines  of  distributed  dislocations  that  are  inclined  to  the 
crack  line.  We  consider  only  the  pile-up  of  edge  dislocations  for  the  mode  I-crack. 
Our  approach  is  quite  different  from  previous  ones  in  that  it  searches  for  the 
possibility  of  existence  of  a  plastic  zone  in  the  form  of  a  linear  array  of  edge 
dislocations  with  a  certain  angle  of  inclination.  In  this  way  the  size  of  the  DFZ 


Introduction. 


and  the  plastic  zone  and  its  position  arrise  naturally  from  the  calculations.  We 
have  also  determined  the  exact  form  of  the  dislocation  distribution.  The  analysis 
is  based  on  a  linear  isotropic  elastic  material. 


Plane  Strain  Solution. 


In  order  to  find  the  force  field  on  an  edge  dislocation  due  to  a  crack  and  other 
defects  it  is  necessary  to  determine  the  stress  field.  For  our  study  the  general  form 
of  the  stress  field,  arising  from  an  arbitrary  symmetrically  loaded  crack  in  plane 
strain,  is  needed.  This  problem  can  be  split  up  into  a  symmetric  mode  I  problem 
and  an  anti-symmetric  mode  II  problem. 

The  Boundary  Conditions  (B.C.)  of  the  symmetric  mode  I  problem  are: 


°xyly  =  0  ”  0 

(1) 

^yylcrack  —  P(*) 

(2) 

^ylouuide  crack  ~  ® 

(3) 

p(x)  =p(-X) 

(4) 

Plane  Strain  Solution. 
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where  a*/,  (k,l  =  x,y)  denote  the  stress  tensor,  (ux ,  Uy)  the  displacement  vector 
#  and  p(x)  is  the  traction  applied  to  the  crack  surface. 


The  stress  field  for  this  problem  is  given  by  (Sneddon  0°>  p  26): 


°w(**y)  =  v/nJ0  (y5-  O4v(4)c_tycos(4x)d^ 

(5) 

« 

Oyy(*.y)  =  -  ^/f  Jo  (y£  +  l)4v(^)e_4ycos(4x)d4 

(6) 

oxy(x,y)=  -yv''|jo  42v(4)e“^ysin(^x)d^ 

(7) 

• 

°zz(x-  y)  =  v(°xx  +  °yy) 

(8) 

• 

where 

v(-)  =  vH'^t)Jo(’t)dt 

(9) 

• 

„  2t  r>  p(u)du 

f(t)  =  - 

n  Jo  1.2  2 

yJX  -  U 

(10) 

• 

Here  Jn(z)  denotes  the  Bessel's  function  of  the  first  kind  of  order  n. 

The  B.C.  of  the  antisymmetric  mode  11  problem  are: 

* 

°yyiy  =  0  ® 

(11) 

• 

°xylcr»ck  —  —  *l(*) 

(12) 

Plane  Strain  Solution. 
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to 


'bJ  outside  crack  0 

(13) 

11 

1 

? 

it 

(14) 

where  q(x)  is  the  shear  stress  applied  to  the  crack  surface. 

The  solution  of  this  problem  is: 

Oxx(x,y)  =  "Vnij  (4y  “  2)£v(£)e~*ycos(!;x)d!; 

(IS) 

Oyy(x.  y)  =  V  n  J0  ^2yv(^e  "  ky  cos(4x)d^ 

(16) 

Oxy(x.y)  =  v''i  Jo  (4y  -  l)^v(4)e  ^y  sin(^x)d^ 

(17) 

07Z(3t,>)  =  VCO^  +  Oyy) 

(18) 

where 

V($)  =  Vf  jof(t)J,(4t)dt 

(19) 

m-l!' uqlu),lu 

"V2-*2 

(20) 

These  solutions  involve  quite  tedious  triple  integrations.  By  expanding  p(x)  and 
q(x)  in  Fourier  series,  we  managed  to  reduce  the  triple  integrations  to  single  ones, 
resulting  in  the  following  solution: 


Plane  Strain  Solution. 
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Symmetric  mode  I: 


00 

P(x)  =  Z  Pn  cos(nnx) 

n  "  0 

or 

cki  =  Z  PnCTki(x- y)  (k,l  =  x,  y) 

n-  0 

where 


o«(x->’)  =  J  *  tJ0(nnt)[yh(x,  y,  t)  -  f(x,y,t)]dt 
o 

<C(*.y)  =  -  J*  tl0(nnt)[f(x,y,  t)  +  yh(x,  y,  t)]dt 
0 

<*?>(*.  y)  =  -  y  Jo'  U0(n/It)j(x,  y,  t)dt 
Antisymmetric  mode  II: 

ao 

q(*)  =  Z  «In  sin(nnx)  +  c^x 

n*  1 

oo 

°ki  “  Z  qnCTki(x-  y)  (k,  1  =  x,  y) 

n  -  0 

where 


o«(x.y)  =  -  ^'gnWfy^U.y.*) -  2f’(x,y.t)]dt 

o£y(x.y)  =  yJo'  gnC^Mx.y.tJdt 
o£y(*.y)  =  J^gnO^ybsC^y-1) _  f.fr.y.t)]* 


Plane  Strain  Solution. 
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gn(t)  =  tJ  i  (nnt)  n  s  1 


BoW  =  1  / 2 


The  term  q^x  in  (26)  is  to  account  for  a  non  zero  shear  stress  at  the  crack  tip. 
The  functions  Ax<y>t)<  Kx<y,  0.  j(.x,y,t),  f(x,y,t),  h*(x,y,t),  fs(x,y,t)  and 
fk(x,y,t)  are  collected  in  the  appendix. 


Along  the  crack  line  a  great  deal  of  simplifications  are  achieved: 


Symmetric  mode  I 


fi  U0(nnt)dt 

«*•»>  -  *1  - 


’  (*’  -  t!) 


OyyCx,  0)  =  o^x,  0) 


0"y(X,  0)  =  0 


Antisymmetric  mode  II 


o"  (x,0)  =  0 


a^(x,0)=  0 


rxny(x,°)=f 


i  tf(t)dt 


°  (x2  -  t2) 


Plane  Strain  Solution. 
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where 


gnW  =  tJ|(nnt)  nil 

(39) 

go(t)  =  t2/2 

(40) 

Hence  we  have  the  solution  of  a  rather  general  problem  in  a  form  attractive  to 
computations.  The  generality  will  allow  us  later  on  to  take  into  account  the 
influence  of  any  kind  of  in-plane  defects  on  the  plastic  zone  size  and  position. 


Method  of  approach. 


Experimental  observations  indicate  that  from  the  crack  tip  emanate  one  or  two 
sets  of  edge  dislocations  making  an  average  angle  ±  a  with  the  crack  line.  At  first 
we  look  at  the  case  of  just  one  set  of  dislocations.  We  idealize  this  picture  by 
assuming  that  the  edge  dislocations  constitute  two  straight  lines  (fig  1). 

We  can  distinguish  four  kinds  of  forces  on  any  of  these  dislocations: 

— * 

1.  The  force  F c  due  to  the  crack  load. 

Fc(t)  =  b(t)  F0(t)  (41) 

Method  of  approach.  9 


where  b(t)  is  the  Burger's  vector  density  at  t, 


F„(t)  =  {[oJy(t)cos(a)  +  ay7(t)sin(a)]'e  . 
-  [o^(t)  cos(a)  +  o*y(t)sin(a)]‘ey} 


°*/(0  (k,l  =  x.y)  is  the  stress  field  at  t  due  to  the  crack  load  and  t  is  a  local 
coordinate  along  the  plastic  zone  (/e[0,Z,]),  fig  1. 


2.  The  force  F  4  due  to  the  other  dislocations. 


?d(t)  =  b(t)jLK1(t,;)b(;)dC 


where 


K  ,(!.?)  = 


2n(l  -  v)(t  -  ;) 


Here  p  is  the  shear  modulus  and  v  is  Poisson's  ratio. 


3.  The  force  F  ctj  due  to  the  interaction  of  the  dislocations  with  the  crack. 


FalW  =  Mt)jLb(;)K2(t,C)dC 


where 


Method  of  approach. 
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K2(t.£)  =  |[o“(t,;)cos(a)  +  o“.(t,  C)  sin(a)J  ex 
+  o^Xt,;)sin(a)]'e  y| 


acifi{t , 0  is  the  stress  field  at  t  due  to  the  interaction  of  the  crack  with  a  unit 
delta- Dirac  dislocation  density  at  £.  C  and  t  are  local  coordinate*  along  the  plastic 
zone  (r,  £  e  [0,  Z.]). 


4.  The  force  F /  due  to  friction  in  the  glide  plane  of  the  dislocation. 


FKO  =  —  Cfb(t)  e  , 


where  Cf  is  of  the  order  of 


cf  =  10  4p  a  10  2p 


The  equilibrium  configuration  of  the  system  requires  that: 


Fc+  Fd  +  Fcd  +  Ff  =  0 


Fo(t)  +  JoL  [k  ,(1,0  +  K2(t,C)]b(Od;  -  cfe  ,  =  0 


At  first  we  neglect  the  interaction  forces.  Hence  we  have: 


?0(t)  +  JLK|(t,;)bQd;-Cfrt  =  o 


/VVV 


Method  of  approach. 
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Due  to  the  chosen  configuration  the  forces  between  the  dislocations  F  4  are  parallel 
to  their  line  of  position,  i.e.  they  have  the  same  angle  of  inclination.  Indeed, 
substituting  (44)  in  (51)  one  obtains: 

F0(t)  +  Tt[joLK,(t,;)b(;)d;  -Cf]  =  o  (52) 

where 

K,(t.O=  |K,(t,C)|  (53) 

Hence,  in  order  to  obtain  equilibrium  the  forces  due  to  the  crack  load  must  also 
have  an  inclination  angle  a  (a  necessary  condition).  From  this  we  deduce  the 
following  approach  to  the  solution  of  our  problem: 

1.  We  choose  an  angle  a. 

2.  We  determine  the  force  field  F  o  due  to  the  crack  load  on  an  edge  dislocation 
with  this  angle  of  inclination. 

3.  In  this  furcv  field  we  look  for  a  straight  line  with  inclination  angle  a,  of  forces 
having  this  same  inclination  angle.  This  gives  us  the  size  and  position  of  the  plastic 
zone. 

4.  Knowing  the  force  field  along  this  line,  we  can  calculate  the  edge  dislocation 
distribution  (Burger's  vector)  which  will  generate  a  force  field  with  the  opposite 
sign  (force  equilibrium).  This  requires  the  solution  of  the  following  integral  equation: 

JLK,(t,0b(0d;  =  -f(t)  (54) 

0 

where 


Method  of  approach. 
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f(t)  =  F 0(t)  -  cf 


F0(t)  =  I  F0(t)| 


5.  The  interaction  forces  between  crack  and  dislocations  can  be  treated  by  iteration. 


Let  the  solution  of  (54)  be  bo{t).  The  dislocation  density  after  one  iteration 


is  given  by  the  solution  of 


?o(t)  +  f  K  ,(t,  ;)b,(4)d;  +  JL  K2(t,  OboM  -  cfe  t  =  0  (57) 


The  iterated  location  of  the  plastic  zone  is  determined  by  the  force  fidd 


F,(t)=  F0(t)  + JLK2(t,;)b0Qd; 


and  the  density  b](t)  is  the  solution  of 


fLK,(t.;)b,(;)d;-  -f,(t) 


where 


fj(t)  =  |F,|  -  cr 


In  case  of  two  symmetrically  located  dislocation  lines  (fig  2)  we  have  to  introduce 


some  additional  forces. 


5.  The  force  F  c<t  due  to  the  interaction  of  the  symmetrically  located  dislocations 


with  the  crack. 


Method  of  approach. 


IV 


* 

I 

I 

i 


F«u(t)  =  b(t)JoLb(;)K3(t,;)d; 


(61) 


where 

K3(t,  0  =  {  [®x>*(ti  0  cos(a)  +  o*s(t,  0  sin(a)]  ~c  x 
-[o«5(t,C)cos(a)  +  o^s(t,;)sin(a)]'ey} 


\ 

i 

t 


(62) 


*>)  is  the  stress  field  at  t  due  to  the  interaction  of  the  crack  with  a  unit 
delta-Dirac  dislocation  density  at  £*.  is  symmetrically  located  with  respect  to  £: 
if  t  =  £  corresponds  to  the  rectangular  coordinate  (x,y)  then  £  corresponds  to  (x,-y). 

♦ 

6.  The  force  F  &  due  to  the  symmetrically  located  dislocations. 


Fd$(t)  =  b(t)JoLb(;)K4(t,Od;  (63) 

where 

K4(t,  C)=|  [oxy(t,  0  cos(a)  +  o*(t,  0  sinfa)]-?  x 

(64) 

“  [a*x(t'Qcos(a)  +  ox!y(t,0sin(a)]‘eyj- 
a$(f,  0  is  the  stress  field  at  t  due  to  a  unit  delta-Dirac  dislocation  density  at  £  . 


Now  equilibrium  requires: 


Fc  + 


Fd+  F 


od 


+  F 


+  F  ds  +  F  d,  =  0 


(65) 


Method  of  approach. 
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We  again  oniy  retain  Fc,  Fa  and  F /.  Iteration  is  now  based  on  the  solution  of 


F0(t)  +  Jql  [ic2(t,o  +  K3(t,;)  +  K4(t,o]b0(C)d£ 


+  JLK1(t,;)b1(Od;+cf'es=0 


Due  to  the  significance  of  the  correction  forces,  iteration  is  more  important  in  the 


symmetric  case. 


The  dislocation  distribution. 


The  determination  of  the  edge  dislocation  distribution  corresponding  to  the  cal¬ 


culated  force  field  amounts  to  the  solution  of  the  following  Cauchy  integral  equa¬ 


tion  for  b(x): 


rL  b(x)ndx 

Jo  2nd  -  v)(*-y)  Uy) 


ft<p(x)dx 

>.  T=J 


where 


The  dislocation  distribution. 


(69) 


<p(x)  = 


jib(x) 


2n(l  -  v) 


The  solution  of  (68)  which  is  zero  at  L  is  given  by  Nlikhlin  (13>  p  130: 


-  f 

n2Jt(L  - 1)  L  0  1  “  ^  0  LL  ) 


1/2 


<P(t) 


V*(L 


(70) 


t  e[0,  L] 


(71) 


The  first  integral  is  a  Cauchy  integral  whose  singularity  can  be  removed  as  follows: 
we  write  the  integral  as  one  centered  about  the  singularity  and  the  rest:  for  f  i 
L/2  wc  get: 


[;(l  -  o] 


1 1 


Jo  t  - 


c 


■fQd;  = 


J*  [(t  -  u)(L  -  t  +  u)]'  2f(t  -  u)^  +  J  ^  _v~ 

- 1  21  IS 


1/2 


(72) 


f(C)dC 


Then  w’c  subtract  the  constant  contribution  of  the  numerator  of  the  first  integral 
which  removes  the  singularity.  In  this  way  we  obtain: 
for  t  £  L/2: 


I  “  J/  {[(t  -  u)(L  -  t  +  u>r«t  -  u)  -  [t(L  -  t)]'  \t)}  ^ 


fL[«L-;)]'g 
>,  .  - 1  ^ 


(73) 


for  t  >  L/2: 


The  dislocation  distribution. 
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(74) 


i  =  Jo  t_c  n;)d; 

+  £"L‘{[(t  -  u)(L  -  t  +  u)]'  2f(t  -  u)  -  [t(L  -  t)],/2f(t)}  ^ 


•  Results  and  discussion. 
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We  ran  the  computer  program  for  a  =  45°,  55°,  60°,  70°  and  80°.  These  angles 
were  chosen  because  experimentally  edge  dislocation  emission  happened  mainly 
in  a  direction  within  45°  to  90°  angles  with  the  crack  plane. 

Figures  3  to  7  give  a  detailed  overview  of  our  scheme  for  one  crack  line  with  a 
=  60°.  Crack  length  and  load  were  taken  unity.  Figure  3  represents  the  force  field 
due  to  the  crack  load  on  an  edge  dislocation  with  angle  60°.  Figure  4  displays  the 
line  along  which  the  forces  are  inclined  over  60°.  The  plastic  zone  is  given  by  that 
portion  of  the  line  that  has  a  line  of  inclination  of  60°(fig  5)  to  within  a  tolerance 
of  4°.  The  force  f(t)  along  the  (normalized)  plastic  zone  is  given  by  fig  6.  Finally 
fig  7  shows  us  the  edge  dislocation  distribution  along  the  plastic  zone.  We  aee 
that  the  method  gives  us  the  plastic  zone  and  DFZ  in  a  natural  way. 

The  effects  of  subsequent  iterations  on  the  plastic  zone  location  are  represented 
in  fig  8.  It  is  clear  that  in  the  case  of  only  one  dislocation  line  iteration  doein  not 
influence  the  results  appreciably. 


Results  and  discussion. 
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Fig  9  points  out  to  the  influence  of  the  friction  on  the  dislocation  density.  As 
long  as  the  friction  force  does  not  exceed  the  minimum  force  density  along  the 
dislocation  line,  the  extent  of  the  plastic  zone  is  not  affected. 


We  also  ran  the  program  for  different  angles,  the  results  of  which  are  displayed 
in  fig  10.  We  notice  that  the  size  of  the  DFZ  is  not  really  dependent  on  the  angle 
of  inclination,  while  the  plastic  zone  size  clearly  is.  There  is  a  general  trend  for 
the  plastic  zone  to  get  shorter  with  greater  angles.  The  forces  along  the  plastic 
zone  and  the  dislocation  distribution  increase  slightly  with  larger  angles. 

The  situation  for  the  two  symmetric  dislocation  zones  is  different.  Fig  1 1  makes 
it  clear  that  iteration  is  quite  important.  The  forces  neglected  in  the  zeroth  iteration 
are  indeed  significant.  Fig  12  displays  an  overview  for  several  angles  of  inclination. 
We  see  that  the  dislocation  zones  are  slightly  smaller  than  the  case  of  a  single 
dislocation  line. 


Summary 


In  this  article  we  studied  the  existence  of  a  plastic  zone  consisting  of  edge  dislo¬ 
cations  around  a  mode  1  cTack.  The  analysis  showed  a  decreasing  plastic  zone  size 
for  an  increasing  angle  of  inclination.  The  presence  of  a  DFZ  emerged  as  a  natural 
product  of  the  calculations.  The  present  method  can  be  applied  to  mode  II 
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problems  and  the  interaction  of  the  crack  with  material  defects.  The  influence  of 
the  generation  of  Helium  bubbles  on  size  and  location  of  the  plastic  zone  will  be 
discussed  in  a  subsequent  paper. 
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Appendix 


We  can  express  J[x,  >-,/).  h{x,y,l),  j{x,y,t),  f  (x,y,l),  h*  (x,y,t),  h,{x,y,  t)  and 
fs(x,y,  i)  as  a  function  of  3  other  functions  H\(p),  Hj(p,  4)  and  H3(x,y,  t ): 


f(*.y.t)  =  ~H,(y) 

(75) 

h(x,  y,t)  =  H2(y,y) 

(76) 

j(x,y,t)  =  H2(x,  y) 

(77) 

f* U.y.t)  =  -  H,(t) 

(78) 
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h  (x,y,t)  =  H2(y,t) 


hs(x,y,t)  =  H2(x,  t) 

f^x.y.O  =  Hj(x,y,t) 

where 

Hi(p) =  ~V(gi  1 2  +^B\l)  r_|gr32giP  +g2Pgi  -  g^gr2 

2v/ 2  L  ^ 

—  1  -  ]  '2  _  |  —  3  2  P  ]  — 3/2 

H2(P-q)  =  ~=(gi  8281  )  [_2g’  ®»p  +g2p8i  -  8ipg2g 

V  2 

[1-32  -  2*1  1  -1/2  —  1  ”1/2 

~2«1  Elq  +g2qgl  ~  glqg2gl  J  +  +  «2»1  )  * 

2g'  ®>Pq  +  (g2pqg)  +  B2pg]q  ~  Elpqg2  —  glpg2q)gl  4" 

^  Elqglp  (g2pgl  glpg2^glqg|  ^ 

P.q  =  x,y  or  t 


H3(x,y,t)  = 


•2x(g3ty  +  &tt  y) 


(g3  +  &»)(g3  +  &a)2  ~  4x21/2 


2*(g3l  +  g*,) 


(g3  +  ga)2((g3  +  &i)2  ~  4*2) 


[(g3y  +  g4y)((g3  +  g*)2  “  4x2)1/2  + 


+  (g3  +  g4)2(g3y  +  g,y)«g3  +  g*)2  -  4x2)_1/2] 


Appendix 


in  which: 


g,  =  (y2  +  t2  -  x2)  +  4x2y2 


2  2  2 
g2  -  y  + 1  -  * 


g3  =  [y2  +  (t  +  x)2] 
&  =  [y2  +  (x  -  i)2] 


!  2 


!  2 


(85) 


(86) 


(87) 


(88) 


The  subscripts  x,  y,  t,  p  and  q  represent  partial  derivatives  with  respect  to  x,  y,  t, 
p,  and  q. 
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